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^ , Abstract 

A mapping is obtained relating radial screened Coulomb systems with 
low screening parameters to radial anharmonic oscillators in iV-dimensional 
space. Using the formalism of supersymmetric quantum mechanics, it is 
shown that exact solutions of these potentials exist when the parameters 
satisfy certain constraints. 
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1 Introduction 

Since the appearance of quantum mechanics there has been continual interest 
in models for which the corresponding Schrodinger equation is exactly solvable. 
Solvable potential problems have played a dual role in physics. First, they rep- 
resented useful aids in modeling realistic physical problems, and second, they 
offered an interesting field of investigation in their own right. Related to this 
latter area, the concept of solvability has changed to some extent in recent years. 
With regards to solvability of the Schrodinger equation there are three interesting 
classes of the potentials. 

The first class is the exactly solvable potentials allowing to obtain in explicit 
form all energy levels and corresponding wave functions. The hydrogen atom 
and harmonic oscillator are the best-known examples of this type. The second 
class is the so-called quasi-exactly solvable (QES) potentials for which a finite 
number of eigenstates of corresponding Hamiltonian can be found exactly in ex- 
plicit form. The first examples of such potentials were given in [1]. Subsequently 
several methods for generating partially solvable potentials were worked out, as 
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a result many QES solvable potentials were found [2]. The third class is the 
conditionally-exactly solvable potentials [3] for which the eigenvalue problem for 
the corresponding Hamiltonian is exactly solvable only when the parameters of 
the potential obey certain conditions. 

Although modern computational facilities permit the construction of solutions 
for any well-behaved potentials to any degree of accuracy, there remains contin- 
uing interest in exact solutions for a wider range of potentials. In connection 
with this, the technique of changing the independent coordinate has always been 
a useful tool in the solution of the Schrodinger equation. For one thing, this 
allows something of a systematic approach, enabling one to recognize the equiv- 
alence of superficially unrelated quantum mechanical problems. For example, 
solvable Natanzon [4] potentials in nonrelativistic quantum mechanics are known 
to group into two disjoint classes depending on whether the Schrodinger equa- 
tion can be reduced to a hypergeometric or a confluent hypergeometric equation. 
All the potentials within each class are connected via point canonical transfor- 
mations. Gangopadhyaya and his co-workers [5] established a connection also 
between the two classes with appropriate limiting procedures and redefinition of 
parameters, thereby inter-relating all known solvable potentials. In order for the 
Schrodinger equation to be mapped into another Schrodinger equation, there are 
severe restrictions on the nature of the coordinate transformation. Coordinate 
transformations which satisfy these restrictions give rise to new solvable problems. 
When the relationship between coordinates is implicit, then the new solution are 
only implicitly determined, while if the relationship is explicit then the newly 
found solvable potentials are also shape invariant [5, 6]. In a more specific special 
application of these ideas, Kostelecky et al. [7] were able to relate, using an ex- 
plicit coordinate transformation, the Coulomb problem in iV-dimensions with the 
iV-dimensional harmonic oscillator. Other explicit applications of the coordinate 
transformation idea can be found in the review articles of Haymaker and Rau [8]. 

Moreover, recently an anti-isospectral transformation called also as duality 
transformation was introduced [9]. This transformation relates the energy levels 
and wave functions of two QES potentials. Many recent papers [10-12, and the 
references therein] have addressed this subject of coordinate transformation plac- 
ing a particular emphasis on QES power-law potentials, which is also the subject 
of the present work in some extent. The generalization to higher dimensions of 
one-dimensional QES potentials was discussed in a recent paper [13] and a few 
specific iV-dimensional solutions were listed. In that work, applying a suitable 
transformation, these potentials were shown to be related to the isotropic oscil- 
lator plus Coulomb potential system and some normalized isolated solutions for 
this system were obtained. 

The importance in the study of QES potentials, apart from intrinsic academic 
interest, rests on the possibility of using their solutions to test the quality of 
numerical methods and in the possible existence of real physical systems that 
they could represent. For instance, anharmonic oscillators and screening Coulomb 
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(or Yukawa) potentials represent simplified models of many situations found in 
different field of physics. These problems have been studied for years and a 
general solution has not yet been found. 

The problem of quantum anharmonic oscillators has been the subject of much 
discussion for decades, both from an analytical and a numerical point of view, 
because of its important applications in quantum- field theory [14], molecular 
physics [15], and solid-state and statistical physics [16, 17]. Various methods 
have been used successfully for the one-dimensional anharmonic oscillators with 
various types of anharmonicities. Relatively less attention has been given to the 
anharmonic oscillators in higher-dimensional space because of the presence of 
angular-momentum states that make the problem more complicated. The recent 
works [10, 11] have shown that there are many interesting features of the an- 
harmonic oscillators and the perturbed Coulomb problems in higher-dimensional 
space, and discussed the explicit dependence of these two potentials. 

Using the spirit of the works in Refs. [10, 11], we show the mappings between 
screened Coulomb potentials (with low screening parameters) and anharmonic 
oscillator potentials in iV-dimensional space, which have not been previously 
linked. The connection between these potentials are also checked numerically by 
the use of the Lagrange-mesh calculation technique [18, 19]. Next we study the N- 
dimensional screened Coulomb problem and higher order anharmonic oscillators 
within the framework of supersymmetric quantum mechanics (SSQM) [20] and 
have shown that SSQM yields exact solutions for a single state only. 

In the following section, we outline the mapping procedure used through the 
article and give the applications. We also discuss in the same section exact super- 
symmetric treatments of the ground state solutions for the initial and transformed 
potentials. Analysis and discussion of the results obtained are given in section 3. 
Section 4 involves concluding remarks. 



2 Mappings between the two distinct systems 



The radial Schrodinger equation for a spherically symmetric potential in N- 
dimensional space 



1 

~2 

is transformed to 
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where ^ (r) = r^ N 1 ^ )//2 i?(r), the reduced radial wave function , and M = N + 2i. 
Note that the solutions for a particular central potential are the same as long as 
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M remains unaltered. For instance, the s-wave eigensolutions and eigenvalues in 
four-dimensional space are identical to the p-wave solutions in two-dimensions. 

We substitute r = ap 2 /2 and R = F(p)/p x , suggested by the known transfor- 
mations between the Coulomb and harmonic oscillator problems [7, 21] and used 
by [10, 11] to show the mappings between unperturbed Coulomb and anharmonic 
oscillator problems, and transform Eq.(l) to another Schrodinger-like equation 
in N' = 2N — 2 — 2A dimensional space with angular momentum L = 2£ + A, 



d 2 F N' - 1 dF 
+ 



dp 2 p dp 



+ LJL^) F=[ ,_ mF (3) 



where 



E - V(p) = Ea 2 p 2 - a 2 p 2 V{ap 2 /2) (4) 

and a is a parameter to be adjusted properly. Thus, by this transformation, 
the N -dimensional radial wave Schrodinger equation with angular momentum I 
can be transformed to a N' = 2N — 2 — 2X dimensional equation with angular 
momentum L = 2£ + A. The significance of the mapping parameter A is discussed 
below. 

A student of introductory quantum mechanics often learns that the Schrodinger 
equation can be solved numerically for all angular momenta for the screened 
Coulomb and anharmonic oscillator problems. Less frequently, the student is 
made aware of the relation between these two problems, which are linked by a 
simple change of the independent variable discussed in detail through this sec- 
tion. Under this transformation, energies and coupling constants trade places, 
and orbital angular momenta are rescaled. Thus, we have shown here that there 
is really only one quantum mechanical problem, not two. 

The static screened Coulomb potential 

Vscij) = -e 2 exp(-5r)/r (5) 

where 5 is a screening parameter, is known to describe adequately the effective 
interaction in many-body environment of a variety of fields such as atomic, nu- 
clear, solid-state and plasma physics. In nuclear physics it goes under the name 
of the Yukawa potential (with e 2 replaced by another coupling constant), and in 
plasma physics it is commonly known as the Debye-Huckel potential. Eq. (5) 
also describes the potential of an impurity in a metal and in a semiconductor [22]. 
Since the Schrodinger equation for such potentials does not admit exact analytic 
solutions, various approximate methods [23, 24, and the references therein], both 
analytic and numerical, have been developed. 

Considering the recent interest in various power-law potentials in the litera- 
ture, we work through the article within the frame of low screening parameter. 
In this case, the screened Coulomb potential appears 

. . , exp(-5r) e 2 2r e 2 5 2 e 2 5 3 2 e 2 5 4 , e 2 5 5 , 
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11 + A 2 + A 3 r + A A r 2 + A 5 r 3 + A 6 r 4 (6) 



r 



in the form. The literature is rich with examples of particular solutions for such 
power-law potentials employed in different fields of physics, for recent applica- 
tions see Refs. [25, 26]. At this stage one may wonder why the series expansion 
is truncated at a lower order. This can be understood as follows. It is widely 
appreciated that convergence is not an important or even desirable property for 
series approximations in physical problems. Specifically, a slowly convergent ap- 
proximation which requires many terms to achieve reasonable accuracy is much 
less valuable then a divergent series which gives accurate answers in a few terms. 
This is clearly the case for the screening Coulomb problem [27]. This also explains 
why leading orders of the 1/N expansion converge at a similar rate for the hydro- 
gen atom, the screening Coulomb potential, and two-electron atom even though 
the last two of these series diverge eventually [28]. In addition, the complexity 
of calculating higher order terms in the series for the corresponding transformed 
potential grows rapidly. Hence, if an accurate approximation cannot be achieved 
in a few terms, the present method may not be useful. However, in section 3 we 
show that the present technique gives excellent estimates for the energy eigenval- 
ues of both, the truncated screening Coulomb and anharmonic oscillator problems 
in higher dimensions. 

Though the mapping procedure introduced is valid for any bound state, 
throughout this work we are concerned ony with the ground state. We have 
two reasons for this restriction. The first reason is that the exact analytical 
ground state solutions for the potentials of interest are available, which provides 
a test for the reliability of the present technique. The second reason is that the 
present model works well for low lying states, which will be shown in section 3. 
Proceeding, therefore, with the choice of a 2 = l/\E n=0 \ in Eq. (4), the screened 
Coulomb problem in Eq. (6) is transformed to an anharmonic oscillator problem 
such that 

T>/ \ / i ^2 \ 2 ^3 4 ^4 6 , ^5 g Aq 1q 
V (P) = 1 + Ttti f P +7^ iTHP +7TF^ r,P+^ IZnP+TT^ ^P 



|£ n =o|r 2|£ n=0 | 3 /^ 4|£ n=0 r 8|£ n=0 | 5 / 2F 16|£ n=0 | 3 
with the eigenvalue 



(7) 



En=o - rj^ (8) 

Thus, the system given by Eq. (6) in iV-dimensional space is reduced to 
another system defined by Eq. (7) in N' = 2N — 2 — 2A dimensional space. 
In other words, by changing the independent variable in the radial Schrodinger 
equation, we have been able to demonstrate a close equivalence between the 
screened Coulomb potential and anharmonic oscillator potentials. 

For almost two decades, the study of higher order anharmonic potentials has 
been desirable to physicists and mathematicians in understanding a few newly 
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discovered phenomena such as structural phase transitions [29], polaron formation 
in solids [30], and the concept of false vacuum in the field theory [31]. Unfor- 
tunately, in these anharmonic potentials, not much work has been carried out 
on the potential like the one defined by (7) except the works in Refs. [32-34]. 
Our investigation in iV-dimensional space, beyond showing an explicit connection 
between two distinct systems involving potentials of type (6) and (7), would also 
be helpful to the literature regarding the solutions of such potentials in arbitrary 
dimensions due to the recent wide interest in the lower-dimensional field theory. 

Eqs. (7) and (8) are the most important expressions in the present work. To 
test explicitly if these results are reliable, one needs to have an exact analytical 
solutions for the potentials of interest, which are found below within the frame 
of supersymmetric quantum mechanics. 

2.1 Supersymmetric treatment for the ground state 

Using the formalism of SSQM [20] we set the superpotential 

W(r) = — + a 2 + a 3 r + a A r 2 , a A < (9) 

for the potential in (6) and identify V + (r) supersymmetric partner potential with 
the corresponding effective potential so that 

V + (r) = W 2 (r) + W\r) = ?^ + [a 2 2 + a 3 (2a 1 + l)}+2(a 1 a 4 + a 4 + a 2 a 3 )r 

r 

+ (2a 2 a 4 + a 2 3 )r 2 + 2a 3 a 4 r 3 + a\r A + Ql ^ ^ 

^ + 2A 2 + 2A 3 r + 2A A r 2 + 2A 5 r* + 2A 6 r 4 ) + (M ~ ~ 3) - 2£ n=0 (10) 

where n — 0, 1,2, ... is the radial quantum number. The relations between the 
parameters in (10) satisfy the supersymmetric definitions 

M — I 2A X A 5 



ai = -J-, « 2 = j^, a 3 = --j=, a 4 = (11) 

Note that in order to retain the well-behaved solution at r — > oo we have chosen 
the negative sign in a 4 . The potential in (6) admits the exact solutions 

*„=o(r) = iV r ai exp [a 2 r + ^r 2 + ^r 3 ) (12) 



where iVo is the normalization constant, with the physically acceptable eigenval- 
ues 

E n=0 = A 2 — - 



4A 2 A 5 

5 -.M 



(M-l) 2 ^/2A & 



(13) 
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under the constraints 



, , , , . 8A 6 A 4 - 2Ai A r— 



(M+l) Ax A £ 



+ 



(M-l)A 



, (14) 



from which one arrives at the uniquely determined values of M ~ 5 and 5 ~ 
0.28 in case of using atomic units in (6). The results obtained here are the 
generalization of the work in Ref . [25] . 

For the anharmonic oscillator potential in (7), we set the corresponding su- 
perpotential 



which leads to 



W(p) = ap 5 + bp 3 + - + dp, a < 0, d < (15) 



*n=o(p) = C oP c exp fo 6 + ^p 4 + d -p>\ (16) 



with Co being the corresponding normalization constant. Note that ^f(p) satisfies 
a differential equation analogous to Eq. (2) and is related to F(p) in Eq. (3) 
as ty(p) = p( N _1 ^ 2 F(p) like the relation between R{r) and ty(r). Identifying 
V + (p) with the effective potential we arrive at an expression for the physically 
meaningful ground state eigenvalues of the anharmonic oscillator potential in 
arbitrary dimensions, 



E - d (2c I n_ 8 A^~2Ag M' 




where M' — N' + 2L, and the relations between the potential parameters satisfy 
the super symmetric constraints 

8 |K =0 | 3 / 2 ' 8a\E n=0 \^' C 2 '° 2a\2\E n= tf °) 

(18) 

As we are dealing with a confined particle system, the negative values for 
a and d would of course be the right choice to ensure the well behaved nature 
of the wave function behaviour at infinity. Our results are in agreement with 
the literature existing for three-dimensions [32, 33] (in case N' = 3) and for 
two- dimensions [34] (in case N' = 2, L — > L — 1/2). 

In sum, we have shown that SSQM yields exact solutions for a single state 
only for the underlying quasi-exactly solvable potentials in higher dimensions with 
some constraints on the coupling constants. These constraints differ from each 
eigenvalue, and hence various solutions do not correspond to the same potential 
and are not orthogonal. We have not found these solutions discussed in the 
literature. 
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2.2 Significance of the mapping parameter 



To show explicitly the physics behind the transformation described above, and 
to make clear the significance of the mapping parameter A, we consider Eqs. (13) 
and (17) together within the same frame and arrive at 

M' A 1 

n= ° = ~M-l\E n=0 \^ (19) 

To be consistent with Eq. (8) we must impose < A < 1 as an integer, such 
that 

M' = 2(iV-l-A) + 2(2l + A) = 
M -I N + 2£-l 1 ' 

It is a general feature of this map that , in both cases (A = 0, 1), the spectrum 
of the iV-dimensional screened Coulomb problem is related to half the spectrum 
of the iV'-dimension anharmonic oscillator for any even integer N'. The reader 
is referred to [7] for a comprehensive discussion of similar conformal mappings in 
the language of physics. 

It is worthwhile at this stage to note that recently Chaudhuri and Mondal [10] 
studied the relations between anharmonic oscillators and perturbed Coulomb po- 
tentials in higher dimensions but their results correspond only to the case when 
A = 1, in this case the three-dimensional perturbed Coulomb problem and the 
four-dimensional anharmonic oscillator cannot be related. However, by intro- 
ducing an extra degree of freedom for the mapping parameter (A = 0) through 
our equations, we can reproduce the well-known results found in the literature 
in three-dimensions. With this exact correspondence we can check Eq. (8), us- 
ing exact results for the screened Coulomb potential, and calculated numerical 
results for the anharmonic oscillator potential. 



3 Results and Discussion 

In this section numerical applications of the transformation presented in the 
previous section are given. Calculations to check the validity of the equations 
developed for the screened Coulomb and anharmonic potentials are also given 
here. Table I displays the exact eigenvalues of the screened Coulomb poten- 
tial in three-, and five- dimensions obtained using the Lagrange- mesh calculation 
technique [18, 19] for selected values of the potential parameters. Highly accurate 
Lagrange-mesh calculation results agree well with the best existing numerical and 
theoretical values obtained in three-dimensions [23, 24]. Due to the constraint in 
the potential parameter A\ expressed in Eq. (14), we are not able to show in the 
same table the corresponding exact energy values which can be calculated by Eq. 
(13). For the work of interest in this paper we set A 1 = —1, consequently the 
adequate S- values satisfying the condition in Eq. (14) fall outside the scope of 
the presented work which has been performed for only low screening parameters. 



8 



Further, our calculation results shown in Table I make clear that the eigenval- 
ues of the five-dimensional screened Coulomb problem with any angular momen- 
tum quantum number £, for a particular 5 value, are equal to the same system 
with £ + 1 in three- dimensions, due to M = N + 2£ which remains unaltered for 
these states. The tabulated results support the work of Imbo and Sukhatme [35] 
in which they formulated SSQM for spherically symmetric potentials in N spatial 
dimensions and showed that the supersymmetric partner of a given potential can 
be effectively treated as being N + 2 (£—>£ + 1) dimensions. This fact was 
exploited in their calculations using the shifted 1/JV expansion. 

It is also noted that for very small values of the screening parameter, the 
screening Coulomb potential reduces to the Coulomb potential that is shape in- 
variant having supersymmetric character. Therefore, the related supersymmetric 
partner potentials, such as Ve and Vi+i, are expected to have the same spectra 
except the ground-state energy. This can easily be seen in Table I for the case 
of 5 = 0.001 in both arbitrary dimensions. For instance, the supersymmetric 
partner of the s-orbital (£ = 0) spectrum of hydrogen is the p-orbital (£ — 1) 
spectrum of the same system. 

Finally, the exact calculated eigenvalues, by the use of Eq. (8), for the anhar- 
monic oscillator in four-dimensions from the known exact results for the screened 
Coulomb problem in three-dimensions are displayed in Table II. These eigenval- 
ues are checked by the Lagrange-mesh calculations and tabulated in the same 
table. For low lying states, the results obtained with the present technique agree 
well with the numerical calculations, but this agreement deteriorates quickly for 
higher lying states. 

4 Concluding Remarks 

The mapping problems in arbitrary dimensions have been the subject of several 
papers and have served to illustrate various aspects of quantum mechanics of 
considerable pedagogical value. As the objective of this presentation we have 
highlighted a different facet of these studies and established a very general con- 
nection between the screened Coulomb and anharmonic oscillator potentials in 
higher dimensional space through the application of a suitable transformation, 
the purpose being the emphasize the pedagogical value residing in this interrela- 
tionship between two of the most practical applications of quantum mechanics. 
The exact ground state solutions for the potentials considered are obtained an- 
alytically within the framework of supersymmetric quantum mechanics, which 
provides a testing ground for benchmark calculations. Although much work had 
been done in the literature on similar problems, an investigation as the one we 
have discussed in this paper was missing to our knowledge. 
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Table 1: The first four eigenvalues of the screening Coulomb potential in Eq. (6) 
as a function of the screening parameter in atomic units. 
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-0.005 891 




2 


-0. 


.037 


515 


-0. 


.014 


939 


-0. 


005 


653 


-0.001 521 


0.025 





-0. 


.475 


461 


-0 


.101 


776 


-0. 


034 


329 


-0.012 495 




1 


-0. 


101 


492 


-0 


.034 


079 


-0. 


.012 


287 


-0.003 770 




2 


-0. 


033 


573 


-0 


.011 


865 


-0. 


.003 


458 


0.000 253 



In five- dimensional space 



0.001 





-0. 


124 


002 


-0 


.054 


561 


-0.030 


261 


-0.019 018 




1 


-0. 


.054 


561 


-0 


.030 


260 


-0.019 


017 


-0.012 914 




2 


-0. 


,030 


259 


-0 


.019 


016 


-0.012 


912 


-0.009 237 


0.005 





-0. 


120 


062 


-0 


.050 


708 


-0.026 


526 


-0.015 428 




1 


-0. 


,050 


684 


-0 


.026 


503 


-0.015 


406 


-0.009 474 




2 


-0. 


026 


468 


-0 


.015 


373 


-0.009 


443 


-0.005 961 


0.010 





-0. 


.115 


245 


-0 


.046 


153 


-0.022 


313 


-0.011 622 




1 


-0. 


.046 


061 


-0. 


.022 


228 


-0.011 


543 


-0.006 070 




2 


-0. 


022 


099 


-0 


.011 


425 


-0.005 


965 


-0.002 980 


0.020 





-0. 


.105 


963 


-0 


.037 


852 


-0.015 


232 


-0.005 891 




1 


-0. 


.037 


515 


-0 


.014 


939 


-0.005 


653 


-0.001 521 




2 


-0. 


.014 


491 


-0 


.005 


286 


-0.001 


263 


0.000 885 


0.025 





-0. 


.101 


492 


-0 


.034 


079 


-0.012 


287 


-0.003 770 




1 


-0. 


033 


573 


-0 


.011 


865 


-0.003 


458 


0.000 253 




2 


-0. 


Oil 


216 


-0 


.002 


974 


0.000 524 


0.003 087 



13 



Table 2: Ground-state eigenvalues of the anharmonic potential in Eq. (7) 



In four-dimensional space 



6 


I 


L 


(taken from Table I) 


Lagrange-mesh calculations 


E n =o 

Exact value (Eq. 8) 


0.001 








0.499 000 


2.831 259 


2.831 259 




1 


2 


0.124 002 


5.679 579 


5.679 573 




2 


4 


0.054 561 


8.562 285 


8.562 268 


0.005 








0.495 019 


2.842 624 


2.842 622 




1 


2 


0.120 062 


5.772 014 


5.772 012 




2 


4 


0.050 684 


8.883 704 


8.883 714 


0.010 








0.490 075 


2.856 927 


2.856 924 




1 


2 


0.115 245 


5.891 401 


5.891 406 




2 


4 


0.046 061 


9.318 882 


9.318 871 


0.020 








0.480 296 


2.885 862 


2.885 862 




1 


2 


0.105 963 


6.144 014 


6.144 024 




2 


4 


0.037 515 


10.325 883 


10.325 891 


0.025 








0.475 461 


2.900 499 


2.900 498 




1 


2 


0.101 492 


6.277 884 


6.277 896 




2 


4 


0.033 573 


10.915 282 


10.915 281 
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